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1.1 Introduction 

The plasma chemistry of reactive gas discharges depends strongly on the collision (elemen- 
tary) processes operating, for given external control parameters, between the constituents of 
the discharge: electrons, ions, molecules, and atoms. Mathematically, these processes are en- 
coded in the collision integrals of the Boltzmann-Poisson equations (see Chapter XY), which 
are classical equations describing the discharge on the macroscopic length scales defined by 
the mean-free paths of the various species and the screening length. The dynamics of the 
elementary processes, however, takes place on the much shorter microscopic scale. It is thus 
controlled by quantum mechanical principles which will be discussed in this chapter. 

From the quantum-mechanical point of view, see Fig. 11.11 for an illustration, collisions 
proceed through temporary compound states (1) which form and fragment on the atomic 
length and time scales Q, in contrast to the macroscopic scales on which typical plasma 
phenomena such as sheath formation or wave propagation occur. That the scales required to 
complete the momentum, energy, and particle Q transfer (rearrangement) defining a collision 
are well separated from the plasma scales has two profound implications: (i) cross sections 
measured, for instance, in crossed beam or in swarm experiments, can be used without 
modifications as input data for plasma modeling and (ii) theoretical calculations of cross 
sections can neglect the plasma environment, which is a tremendous simplification. 

^If not stated otherwise, in this chapter it is implicitly assumed that physical quantities are measured in atomic units, that 
is, lengths, masses, and charges are given in terms of the Bohr radius as, the electron mass rrte, and the elementary charge e, 
respectively. The unit of energy is then 2Rq, where Rq is the Rydberg energy. 

■^The term "particle " denotes either a single "fundamental" particle, having no internal degrees of freedom relevant for the 
collision process, or a collection of fundamental particles in a bound state; in the latter event the term "fragment" will be also 
used. 



1 



2 



CHAPTER 1. ELEMENTARY PROCESSES IN GAS DISCHARGES 



A +B 




electrode 2 



Figure 1.1: Schematic representation of the length scales in the bulk of a typical low- 
temperature gas discharge {ue = rii = 10^° cm~^, Tg = 10 eV, Ti = 0.03 eV, and 
^gas = 10^^ cm~^). The macroscopic scales on which the Boltzmann-Poisson system op- 
erates are the spatial extension of the plasma L ~ 0{cm), the mean-free-path of the species 
/ ~ O{10~^ cm), and the plasma screening length A ~ O{10~^ cm). On the rhs is an 
illustration of how the (temporary) collision compound AB~ , whose length scale is a few 
atomic units, that is, of the order of 10~^cm, controls the first six reactions of Table 11.11 
The branching of the AB~ state at a given energy determines the cross sections for the 
reactions. The double arrow indicates that the fragments can appear in the entrance and 
the exit channel of a collision. 

Whereas the physics of a collision is unaffected by the plasma, the reverse is of course not 
true. Elementary processes determine to a large extend the properties of a plasma (see previ- 
ous chapters). In the first place, a gas discharge can be only maintained electrically, because 
electron impact ionization and, for an electro-negative gas, dissociative electron attachment 
produce positively and negatively charged carriers. Equally important are electron impact 
excitation and dissociation of molecules. They transfer not only external electric energy into 
internal energy but, most importantly, they also produce the species which are eventually 
utilized in the technological application of the discharge: electronically excited species when 
the discharge is used as a light source or laser and reactive fragments (radicals) when the 
discharge is employed for surface processing or catalysis. 

Technologically interesting gas discharges contain complex molecular gases such as CF4, 
CF3I, CsFs, CCI2F2, or SFq with a multitude of excited and fragmented species. Even 
simple diatomic molecules, for instance, O2 or N2 give rise to a large number of species, with 
an accordingly large number of elementary processes. Leaving aside elementary processes 
containing three particles in the entrance channel, which are only relevant at rather high 
densities, the most common collisions for a generic electro-negative, diatomic gas are shown 
in Table 11.11 where they are also classified according to the collision compounds controlling 
the microphysics: AB~ , the compound illustrated on the rhs of Fig. II. H A"^" , A2B, A2B~ , 
and AB. The branching of the compounds and thus the probabilities for the various collisions 
(cross sections), depend on the initial energy and on the properties of the compound states 
at the distance where they have to lock-in into the asymptotic scattering states defining the 
various collisions. 

In view of the great technological and economical impact of reactive gas discharges, it 
is somewhat surprising, that the number of experimental groups measuring cross sections 
for gases of plasma-chemical relevance is rapidly diminishing. Thus, any listing of currently 
available cross section data must be necessarily incomplete. Since, in addition, the gases 
of interest change with time, cross section data are not included at all in this chapter. As 
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Table 1.1: Typical collisions and their compound states for an electro-negative, diatomic 
gas. The indices i and / denote initial and final internal states of the molecules, atoms, and 
ions before and after the collision, respectively. 



collisions with 

(1) e- + AB 

(2) e" + AB 

(3) e" + AB 

(4) e" + AB 

(5) A- + Bi ^ 

(6) A- + Bi ^ 
collisions with 

(7) e- + Ar - 

(8) Ar + AB+ 

(9) Ar + AB+ 

(10) Af + ABi 

(11) Af + ABi 

(12) e- + AB+ 



compound state AB^ 



> AB- 

> AB- 

> AB- 

> AB- 
AB- 
AB- 



A 



e" + ABf 
2e- + AB+ 
f Bf 

- Af + Bf 
Af + Bf 
ABf 



different compound states 

A2- ^ 2e- + Af 
A2B ^ Af + ABf 
A2B ^ Af + Af + Bf 
A2B" -> e~ + Af + ABf 
A2B- ^e- + AsBf 

^ AB ^ Af + Bf 



elastic and inelastic scattering 
electron impact ionization 
dissociative attachment 
dissociation 
direct detachment 
associative detachment 

impact detachment 
ion-ion annihilation 
ion-ion annihilation 
direct detachment 
associative detachment 
dissociative recombination 



far as they exist, they can be found, for instance, in the review article by Brunger and 
Buckman (j2|), the monograph by Christophorou and Olthoff ([3[), and in web-based cross 
section compilations sponsored by national research institutions and university groups whose 
activities depend on atomic and molecular collision data. The largest ones, maintained, 
respectively, by the International Atomic Energy Agency, the American National Institute of 
Standards and Technology, by Oak Ridge National Laboratory, USA, the Japanese National 
Institute for Fusion Science, the Weizmann Institute of Science in Israel, and the Universite 
Paris-Sud in Orsay, France, are (j3): 



1. http://www.iaea.org/programmes/amdis 



2. http: / / physics.nist.gov/PhysRefData/contents.html 



3. http: / / www-cfadc.phy.ornl.gov 



4 . http : / / dbshino .nifs.ac.jp 



5 . http : / / plasma-gate . weizmann .ac.il/ 



6. http: / / gaphyor.lpgp.u-psud.fr 



Due to lack of empirical cross section data, theoretical calculations become increasingly 
important. There are two major lines of attack. One, the ab-initio approach, either attempts 
to construct parameter-free eigenstates for collision compounds and to match these states 
to the asymptotic scattering states representing specific collision products or to directly 
obtain parameter-free scattering matrices from a suitable variational principle ([6|). Both 
approaches are extremely time consuming and, so far, have been applied only to a few 
selected collisions involving molecules with at most three atoms. Even when the formal 
problems preventing the application to complex molecules are solved in the future, it is hard 
to imagine a non-expert routinely using an ab-initio package to generate cross sections for 
as yet unstudied collision processes. 
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The other, perhaps more valuable approach, as far as plasma-chemical applications are 
concerned, is the semi-empirical app roach. Examples of which are the binary-encounter (0; 
H; H) and the Deutsch-Mark (ITol : llll: fl^ model for electron impact ionization, the resonance 
model for dissociative attachment (flSl : [T3), electron-detachment ( ITsI : ITgI ). dissociative re- 
combination (jlTl : [iSi), and vibrational excitation (114 : [l9l ). and the Landau-Zener model for 



ion-ion annihilation (l20l : |21| ). Common to all these models is that they attempt to encap- 



sulate the complicated collision dynamics in a few parameters which can be either obtained 
from electronic structure calculations or from experiments. To some extend these models 
are ad-hoc, but they have the great virtue to provide clear physical pictures for the collision 
processes, thereby helping the plasma physicist to develop an intuitive understanding of the 
microphysics of the discharge. 

Despite the great diversity of elementary processes in a gas discharge, this chapter tries to 
give an uniform presentation of the subject within the framework of multi-channel scattering 
theory 



23 



1241 ) focusing, in particular, on inelastic and reactive collisions. On purpose, 
cross section formulae are not developed to a point where they could be apphed to particu- 
lar collision processes, because this involves mathematically rather technical (approximate) 
solutions of the given equations which moreover have to be worked out case- by-case. The 
details can be found in the original literature. Additional information, in particular with 
regards to formal aspects of calculating cross sections, can be extracted from the review arti- 
cles by Rudge Bardsleyiand Mandl ((SI), Inokuti ([23), Lane dSI), Delos dH), Chutjian 
and coworkers f,3Q.) . Hahn (j3ll ). and Florescu- Mitchell and Mitchellp2i). 



1.2 Fundamental concepts 
1.2.1 Collision cross section 

The kinetic description of a gas discharge is based on a set of Boltzmann equations for the 
distribution functions /^(r, v,t) of the participating species i = 1,2,...,N and those parts 
of Maxwell's equations which are necessary to describe the electro-magnetic driving of the 
discharge. For a capacitively coupled radio-frequency discharge, for instance, the equations 
are 

[5i + vV.-^V.<l'-V.]/.(r,v,t) = J] /f [{/,}], (1.1) 

1 ^ r 

= Vg, / dvfi , (1.2) 

where /f [{/,}] is the collision integral due to process p appearing in the Boltzmann equation 
for species i and depending on the subset {/, } of the distribution functions; $ is the electric 
potential 0. The external driving is in this case simply encoded in the boundary conditions 
for the Poisson equation; and are the charge and mass of species i, respectively, and 
eo is the dielectric constant of the plasma. 

The total number N of species required for the kinetic modeling of a generic electro- 
negative, diatomic gas discharge, can be deduced from Table II. 1[ Clearly, besides electrons 
and molecules of the feed-stock gas AB, atoms A and B, and ions AB~^ and A~ have to be 

•^In this subsection atomic units are not used. 
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taken into account. In general, it is also necessary to include for the molecular and atomic 
constituents some excited (meta-stable) states. Thus, in total, = Nab + Na + Nb + 
Nab+ + Na- + 1 species have to considered, where N^ is the number of excited states kept 
for the atomic or molecular constituent a. 

Each reaction shown in Table 11.11 gives rise to collision integrals in the Boltzmann equa- 
tions of the respective reaction educts (species in the entrance channel) and reaction products 
(species in the exit channel). For elastic scattering, the structure of the collision integral is 



similar to the one originally derived by Boltzmann (|33l). The majority of collisions, how- 
ever, is inelastic and reactive, some involve even more than two reaction products. The 
corresponding collision integrals have then to be constructed from scratch, using elementary 
statistical considerations. 

Take, for instance, electron impact ionization, where a primary electron scatters off a 
molecule and produces a secondary electron and a positively charged ion. Assuming the 
target molecule to be at rest, which is a very good approximation because the molecule is 
much heavier then the electron, the collision integral in the electron Boltzmann equation 



due to this process reads (j3J) 



Ii[fe] = CabI^I dYadw"W\Ya,Yy')f,{Wa) + j dWadY'W\way^v)f,{Ya) 

- I rfvWiy^(v;v',v")/e(v)| , (1.3) 

where momentum conservation is taken into account, the variables r and t in the distribu- 
tion functions are suppressed, and cab is the concentration of the molecules. The function 
W^{'Va',^',^") is the probability for an incident electron with velocity to produce two 
departing electrons with velocities v' and v", respectively. Using energy conservation, it can 
be related to the differential cross section for electron impact ionization 



W (V„; V , V )C/V C/V = -Va^ ^ Q iVa,V ,Q^,QJ 

X 6{Vv'^ + v"^ - u)dQ.'jn'^dv'dv" , (1.4) 



with Ei the ionization energy of the molecule, the electron mass, u = \/v\ — 2Ei/me, 
and q\va,v'\Q!a,Q!^) the differential ionization cross section defined by the relation 

da^ = q\va, v\ Q'^, Q'^)dQ'JQ'^dv" (1.5) 

with < v" < \/vl — 2Ei/me (cf. Fig. 11.21 for notational details). 

The sub-nm physics of impact ionization, that is, the momentum, energy and elec- 
tron transfer during this particular electron-molecule collision, is concealed in the function 
{va,v" This statement holds for all processes. The information required about 
elementary processes reduces thus to a set of (differential) cross sections. The direct mea- 
surement of which is tedious, expensive, and, when meta-stable states are involved, which 
cannot be prepared outside the plasma, sometimes even impossible. It is at this point, where 
theoretical calculations of cross sections - ab-initio or otherwise - have a great impact. 

A systematic solution of the electron Boltzmann equation expands the electron distri- 
bution function in terms of velocity space spherical harmonics. The angles fl'^ and fl" can 
then be integrated out leading to an hierarchy of Boltzmann equations for the expansion 
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Figure 1.2: Velocity space kinematics of electron impact ionization. Upon impact of the 
incident electron with velocity vector two electrons are ejected with velocity vectors v' 
and v", respectively. The solid angles dfl'^ and dQ" corresponding to these two velocities are 
defined with respect to the direction of v^. Because of the large mass difference of electrons 
and molecules, the molecule can be assumed to be at rest. 



coefficients (|34). When the anisotropy of the electron distribution is negligible, the lowest 
order equation suffices. The collision integral of interest is then 



:i.6) 



:i.7) 



which is usually sufficient for a particle-based simulation of the ionization process (see Chap- 
ter XY and Ref. (ET ' 



One more integration would lead to the total ionization cross section, 

/(^a) = / dv"q\Va,v") , 



1.2.2 Formal scattering theory 

Collisions affecting the charge balance and chemistry of gas discharges change the type, 
and sometimes even the number, of scattering fragments. An example is electron impact 
ionization, the process used in the previous section to introduce the concept of a collision 
cross section. It is an inelastic, break-up collision with an electron and a molecule in the 
entrance and two electrons and a positive ion in the exit channel. In comparison to the 
entrance channel, the kinetic energy of the relative motion of the fragments in the exit 
channel is moreover reduced by the ionization energy of the molecule. Thus, a microscopic 
descripti on o f electron-impact ionization, and likewise of many of the other processes listed 
in Table [TTTl cannot be based on simple potential scattering theory (elastic scattering). A 
generalized scattering theory is rather required, capable to account for changes in the internal 
energy (inelasticity), for rearrangement, and for break-up of the scattering fragments. 

The aj 
theory 



ropriate theoretical framework is quantum-mechanical multichannel scattering 



23l : [24 ). To introduce its essential ingredients, two colliding fragments are con- 



sidered. In the center-of-mass frame, the total Hamiltonian of the system is 



H = Trel + Hir,t + V = Ho + V 
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in-coming plane wave 




out-going spherical wave 

Figure 1.3: Schematic representation of the incoming pla ne wa ve and the diverging spherical 
wave as described by the Lippmann-Schwinger equation (ll.9p . 

where T^e\ is the kinetic energy of the relative motion, Hi^t controls the internal degrees of 
freedom of both fragments, and V is the interaction energy between the two. The Lippmann- 
Schwinger equation for the scattering state with the boundary conditions shown in Fig. 11.31 
reads in Dirac's bra-ket notation 







;i.9) 



where the first term denotes the incoming plane wave in the entrance channel. The channel 
state, |$ka) = |0a)|k), satisfies {E — ifo)|^ka)- Thus, E = E^a = + uja, where k is 

the relative momentum, Ua is the internal energy, and fi is the reduced mass of the fragments 
in the entrance channel. 

Quite generally, the scattering amplitude, which in turn determines the differential colli- 
sion cross secti on, i s defined as the amplitude of the out-going spherical wave emerging from 
the rhs of Eq. fll.9p for large inter-particle distances. Hence, in order to find the scattering 
amplitude, Eq. fll.9l) has to be expressed in coordinate representation, which is here specified 
by r, the inter-particle distance, and p, the internal coordinates of both particles, and then 
the limit r — > oo has to be taken. Normalizing continuum states on the momentum scale 
leads then to dH; ET 



13 
/3 



exp\ik ■ r]6a/3 - /(k'/3, ka 



exp [ik'r] 



exp[zk ■ Y\5a(3 + ipl,, 



7,(+) 



0/3 (p) for r — > cxD 



:i.io) 



:i.iii 



where k' is the relative momentum after the collision, and the prefactor in front of the 
diverging spherical wave. 



/(k'/3,ka) = -^($k',|l^|vl/W) 



'1.12) 



is the scattering amplitude. 
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The differential collision cross section is microscopically defined as the ratio of the out- 
going particle flux into the solid angle dQ' originating from the diverging spherical wave 

'^^Ikai^) introduced in Eq. (II. lip and the incoming current density due to the plane wave. 
Hence, 

da = , , , 1.13 

k/ fi 

where is the unit vector in the direction of r. Inserting ''/'k'73 k«('") ^'^^^ quantum- 
mechanical expression for the particle flux, j = {l/{2i^))[ip*'Vip — V'V^/'*], yields jout = 
fc'|/(k'/3,ka)pe,/(^r2) and thus 

dao^^p = j\f{k'P,ka)\'dn' = J^||($k,^|y|M/W)|W , (1.14) 
or, when continuum states are normalized on the energy scale and the identity V^^^^ = t$ 



is used (P;l23), 



da^^^ = , (1.15) 

where t = V + V[E — Hq + ir]]^^t is the transition operator (T- matrix). 

In the derivation of Eq. fll.l2p it was implicitly assumed that r is the inter-particle distance 
in both the entrance and the exit channel. The particles remain therefore intact in the 
course of the c ollisio n. Th ey may only change their internal state. Hence, the cross section 
formulae fll.l4p and fll.lSp can be only applied to elastic (cu/j = Ua) and inelastic {uis ^ uja) 
scattering. 

To obtain the cross section for reactive scattering, the fact has to be included that the 
type of the particles, and hence the relative and internal coordinates, change during the 
collision. Thus, writing the Hamiltonian in the form (11.81) is only adequate in the entrance 
channel. In the exit channel, it is more appropriate to partition the Hamiltonian according 
to the reaction products and write 

H = TU^H[,,^V' = %^V\ (1.16) 

where Tj.g[, -ff/^t' ^^"^ relative kinetic energy, the internal energy, and the interaction 

energy in the exit channel. Then, in addition to Eq. (11.91) . the scattering state I^I'ka'') obeys 
also an homogeneous Lippmann-Schwinger equation Q, 

l^i+h = ^ V'\^^^'^) , (1.17) 

from which the scattering amplitude in the exit channel (reaction amplitude) can be deduced 
by the same procedure as before, except that now the coordinate representation with respect 
to r' and p', the relative and internal coordinates in the exit channel, has to be chosen. In 
this representation, the scattering state becomes a diverging spherical wave for r' — > cxd. The 
prefactor of which (continuum states normalized on the momentum scale), 

/'(k'/?,k«) = -^(<i.;,,,|y'|vi/W^ , (1.18) 



The equation is homogeneous because the incoming wave belongs to a different Hilbert space 
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with /i' the reduced mass in the exit channel and 1$^/^) an eigenstate of Hq, can be identified 
with the reaction amphtude. Hence, the differential cross section for reactive scattering is 
given by 

da.^, = ^j\f'{k'(3,ka)\'dQ' = , (1-19) 

where energy conservation enforces now E = k'^/(2fi') + cj^ = k'^/{2fi) + uj^. Obviously, 
Eq. fll.l9p reduces to Eq. (11.141) for V = V which implies Hq = Hq and thus r' = r, fi' = fi, 
and $' = $. 

If the interaction V in the entrance channel is simpler than the interaction V in the 
exit channel, it may be more convenient to use the adjoint scattering state, (^l/^"/]!) which 
describes an incoming wave in the exit channel. The reaction cross section can then be 
written as 

da.^, = (^^l(^Sl^l<^k.)|W , (1.20) 

which contains V instead of V. 

The formalism described so far is only applicable to binary collisions, that is, collisions 
containing two particles in the entrance and exit channel, respectively. The theoretical 
description of collisions involving three (or more) reaction products (break-up collisions) can 
be also based on a set of Lippmann-Schwinger equations, but the increased dimensionality of 
the relative motion requires a substantial extension of the formalism (see Ref. (|22|)). beyond 
the scope of this introductory presentation. Additional complications arise for identical 
particles in a given channel, or when the interaction is long-ranged, as it is, for instance, the 
case for electron-impact ionization. Ionization cross sections are therefore hardly obtained 
from rigorous calculations. They are usually estimated from less ambitious, semi-empirical 
models to be described below. 



1.2.3 Adiabatic approximation 

The formalism described in the previous subsection will be now applied to electron-molecule 
scattering. To be specific, a diatomic molecule with N electrons is considered. The scattering 
wavefunctions, \l/*^^''(r', r, R), from which the scattering amplitudes are obtained, depend 
then on the coordinates of the target electrons, which are collectively denoted by r', the 
coordinate of the projectile electron r, and the inter-nuclear distance R. In the interest of 
clarity, spin degrees of freedom are suppressed. 

The Lippmann-Schwinger equation for electron-molecule scattering cannot be solved ex- 
actly. Approximation schemes, usually based on an expansion of the scattering state in an 
appropriate basis, are therefore necessary. When the target states, that is, the eigenfunc- 

tions of H^^^ = Tfj + where Tr is the kinetic energy of the nuclei and -ffgf^^ is the 

Hamiltonian for the A^ target electrons at fixed R, are used, the expansion leads to the 
close-coupling approximation. Identifying p with (r',R), the close-coupling approximation 
is basically identical to the channel state representation described in the previous subsection. 
However, this brute-force approach is not well suited. The number of channels (electronic, 
vibrational, and rotational), which are all coupled by the Lippmann-Schwinger equation, is 
too large. In addition, due to the differences in the electronic and nuclear energy scales, the 
rate of convergence is rather pure. 
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intemuclear distance R 

Figure 1.4: Potential energy surface of the initial and final state of electron-impact excitation 
of a diatomic molecule and the relationship of the various energies occurring in the adiabatic 
approximation for the scattering cross section. The kinetic energy of the electron before and 
after the collision is denoted by Cj = A;^/2 and e/ = A;^/2, respectively. 



1.2. FUNDAMENTAL CONCEPTS 



11 



An approximation scheme accounting for the difference in the energy scales, and thus more 
appropriate for electron-molecule scattering, is the adiabatic approximation. It expands the 

scattering state in terms of eigenfunctions of H^^^^^\ that is, eigenf unctions of the fixed- 
nuclei Hamiltonian for the A^-|- 1 electrons {N target electrons and one projectile electron) of 
the compound state AB~, which depend only parametrically on the inter- nuclear distance. 
Since the dynamic coupling due to the perturbation Tr is usually negligible, only one term 
contributes to the expansion. The computational costs are thus substantially reduced. In 
addition, provided the conditions for the validity of the adiabatic approximation are fulfilled, 
the target can be also described adiabatically. Thus, the target states are Born-Oppenheimer 
states, which are of course also much easier to determine than the target states needed in a 
close-coupling calculation. 

The adiabatic approximation is valid far away from excitation thresholds and when the 
collision time is much shorter than the period of nuclear vibration and rotation. It is thus 
mostly applied to non-resonant electron-molecule scattering. As shown by Shugard and 
Hazi (|36i) , the differential cross section for electron impact excitation of a diatomic molecule 
reads for continuum states normalized on the energy scale. 



dCmViJi^rifUfJf — , \TnfVjJf,niVi.JiiS^fi^i)\ dVtf 



:i.2ii 



with Vti and f2/ the solid angles of kj and kj, the electron momenta before and after the 
collision, respectively. The total scattering amplitude. 



:i.22) 



is an average of the fixed-nuclei scattering amphtude, 

,1/2 



^.„^(fi^,r]„R) = I rfrrfr'exp[zk;-r]<|.:/r';R)l^(r,r';R)vl/^-3„^(r,r';R) , (1.23) 



(277)3/2 

over the ro- vibrational states FnujijV) satisfying 



Tr + en{R) 



Fw(R) = 



:i.24) 



where uJnvj is the ro- vibrational energy and en{R) is the potential-energy-surface defined as 
the eigenvalue of the fixed-nuclei Schrodinger equation for the target electrons 



$„(r',R) 







:i.25) 



In view of the discussion of the previous subsection, the structure of the formulae f ll.2ip - 
(11.231) should be clear. The physical meaning of the various energies can be deduced from 
Fig. 11.41 Note, in particular, that the energy £ = + = e„. + -E — oj^ViJi appearing 
in Eq. (11.231) is only the total electronic energy, for the projectile electron and the N target 
electrons, in contrast to the total energy E which includes electronic and nuclear degrees of 
freedom. Ro- vibrational states are characterized by the quantum number n of the electronic 
state in which the nuclear motion takes place, the vibrational quantum number i/, and the 
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rotational quantum number J. Energy conservation implies E — kf /2 + cOmuiJi = k}/2 + 

Within the adiabatic approximation, it is necessary to determine the fixed-nuclei scatter- 
ing amplitude and the ro-vibrational states of the molecule, which in turn depend on the 
potential energy surface of the target molecule. The former can be obtained, as a function of 
R, from the asymptotics of the + 1 electron, fixed-nuclei Lippmann-Schwinger equation, 

which determines \l/^'^j^^.(r, r'; R), while the latter requires, again as a function of R, the 

solution of the N electron problem (11.251) . In both cases, anti-symmetrized wavefunctions 
have to be used because of the indistinguishability of electrons. Thus, even when the nuclear 
motion is split off, the calculation of cross sections for electron-molecule scattering remains 
a formidable many-body problem (l28l ). 



1.2.4 Resonant scattering 

At electron energies of a few electron volts or less the collision time is long compared to 
the period of the inter-nuclear motion and the adiabatic approximation fails. The pro- 
jectile electron is then so slow that it is captured by the molecule giving rise to a bound 
state of the negatively charged molecular ion, which is the collision compound for electron- 
molecule scattering 0. This state interacts with the electron-molecule scattering continuum, 
acquires therefore a finite hfetime, and turns into a quasi-bound (auto-detaching) state. 
Auto-detaching states play a central role in (vibrational) excitation, attachment, recombi- 
nation, and detachment collisions. From a theoretical point of view, all these processes can 
be hence analyzed within a model describing a discrete state (resonance) embedded in a 
continuum of scattering state. 

A particularly elegant derivation of the resonance model for an electron colliding with a 
diatomic molecule containing N electrons has been given by Domcke ( (37l ) who uses many- 
particle Green functions to reduce the -|- 1 electron scattering problem to an effective 
single-electron problem. The reduction is achieved by two projections: First, electronic 
states which are not accessible at the energy considered, are eliminated by introducing an 
optical potential for the incoming electron. Then, in a second step, the fixed-nuclei T-matrix 
is split into a rapidly varying part due to the quasi-bound resonant state and a slowly varying 
background term. The splitting of the T-matrix can be shown to be equivalent to an effective 
single-electron Hamiltonian at fixed-nuclei describing a resonance coupled to a continuum 
of states. At the end, the kinetic energy of the nuclei is included and the electronic degrees 
of freedom are integrated out to obtain an effective Lippmann-Schwinger equation for the 
nuclear dynamics which, with appropriate boundary conditions, can then be used to calculate 
the collision cross sections of interest. 

It is essential for the formalism that the optical potential supports a resonance and that 
the resonance can be extracted from the single-electron continuum such that the scattering 
background contains no spurious resonances. The many electron problem is then com- 
pletely buried in an optical potential, which can be calculated separately using, for instance, 
many-body perturbation theory (l38l ). In principle, the formalism can handle more than one 
resonance as well as electronically inelastic processes (M)- So far, however, it has been only 
applied to electronically elastic collisions involving a single resonance and a single potential 
energy surface for the nuclear motion of the target. 



^The electronic and vibrational properties of negative ions play therefore an important role in low-temperature gas discharges 
with kgTf, ~ 0(1 — lOeV), even when the gas is not electro-negative, that is, when the negative ion is unstable on the plasma 
time scale and therefore irrelevant at the kinetic level. 
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The effective Hamiltonian, visualized in Fig. II. 5[ reads 
H = Ho + V 



Ho 
V 



i) 



TR + Vd 



'd\ + / kdkd^l 



Tr + Vo + 



kdkdQi 



|0d)Vdk(0L+^|+h.c. 



:i.26) 



where is the resonant state and l^j;^'') are the scattering states orthogonal to \4>d)- The 
coupling between the two types of states is given by V^k = {4>d\Hei\4>l^^), where H^i = 
— V^/2 + Vopt is the Hamiltonian for a single electron in the optical potential Vopt] Vd = 
{(f)d\Hei\4>d) + Vb is an operator specifying the potential of the resonant state and Vq is the 
corresponding operator for the potential of the target molecule. 



The Hamiltonian (ll.26p . and thus the Lippmann-Schwinger equation associated with it, 
operate in the combined Hilbert space of the scattered electron and the two nuclei. Pro- 
jecting out the electron, switching to the coordinate representation with respect to the 
inter-nuclear distance, and using out-going wave boundary conditions yields, when rotations 
of the molecule are ignored, a Lippmann-Schwinger equation. 



E + 



2^ZR2 



Vd{R) 



^I^e\R) - J dR'F^^\R,R',E)¥;J{R') = J{R) , (1.27) 



with a kernel 



FW(i?, R', E) = j kdkdnkVdk{R)G^^\R, R', E - e/2)V*^{R') , (1.28) 



where 



G^+\R,R',E) = {R\[E + 



1 d^ 
2ildR^ 



-VoiR) + iv] '\R') 



[1.29) 



is the Green function for the nuclear motion on the potential energy surface of the target, 
Vo(-R) = (-R|Vo|-R), and /i is the reduced mass of the nuclei; Vd{R) = {R\Vd\R) = ed{R) + 
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Vo(-R) is the potential energy surface of the resonant state and Vdk{R) = {R\Vdk\R). The 
inhomogeneity J{R) on the rhs of Eq. (11.271) depends on the boundary conditions. It will be 
discussed below for particular collision processes. 

Equation (11.271) is an effective Lippmann-Schwinger equation for the nuclear dynamics 
in the energy- dependent, nonlocal, and complex potential of the resonant state. To make 

this more explicit, G'^\r,R',E) is expressed in terms of a complete set of target nuclear 
wavefunctions Xv{R)- Employing Dirac's identity and assuming that the interaction between 
the resonance and the scattering continuum depends only on k (isotropic interaction), leads 
to 

(i?, E'- E) = A{R, R'; E) - -T{R, R'; E) (1.30) 

with 

A(fi,ff;E) =4,X: p/<iE'lM^i|(^^4§l^ (1.31) 

and 

T{R,R';E) = STxY,VdE-.AR)xl{R)x.{R!)V:E-.S^) , (1-32) 

V 

where the symbol "P" denotes the principal value of the integral and lOy stands for the 
vibrational energies of the target molecule. 

The inverse of r(i?, R!\ E) is the lifetime of the resonance. As expected, the auto-detaching 
property of the resonance arises from its coupling to the scattering continuum. In principle, 
the resonance has a finite lifetime even in the absence of this coupling because the optical 
potential, Kpt, is complex. Its imaginary part induces therefore a lifetime. This contribution, 
however, is much smaller than the one due to VdE- Hence, it is usually neglected. 

The non-locality of the potential complicates the numerical solution of Eq. (I1.27p . In the 



early applications of the resonance model (jl3l : fl4 : flSl : [l6l : [TTI : [TsI : 1261 ). the nonlocal potential 
was th us replaced by a local one. The local approximation can be obtained from Eqs. (I1.3ip - 
(11.321) by identifying the energy available for the scattered electron with an effective resonance 
energy: E — uj„ ^ Ej.es{R) (jlTi). The completeness of the vibrational target states can then 
be used to obtain A{R, R'; E) = Al{R)6{R - R') and T{R, R'; E) = Tl{R)S{R - R') with 

T,{R) = 87r\VdE...iRm\^ ^l{R) = 47r p| '^^(l^l^ij' ' ^^"^^^ 
which reduces the Lippmann-Schwinger equation to an ordinary differential equation: 

¥,^JiR) = JiR) . (1.34) 



1 . t. 



The scattering cross sections derived from Eq. (I1.27P or Eq. (11.340 are only as good as 
the optical potential, Vopt, the potential energy surfaces Vo(-R) and Vd{R), and the coupling 
function Vdk{R). These quantities have to be obtained from separate calculations, ideally 
using ab-initio techniques of molecular structure theory, or directly from experiments. But 
even then, the cross sections form the resonance model are semi-empirical in the sense that 
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a priori assumptions about the relevance or irrelevance of molecular ion and target states 
have to be made from the very beginning. This weakness of the model, however, leads at 
the same time to its strength: Technically tractable equations, with an intuitive physical 
interpretation, which, in the local approximation, can be even solved analytically with semi- 
classical techniques (i4dl : l4ll : 



1.3 Typical processes 

Now representative elementary processes are discussed in more detail, focusing, in particular, 
on reactive and inelastic collisions, which change the internal energy and composition of the 
scattering fragments. Elastic scattering between the various species is not explicitly included 
although the associated cross sections are usually much larger than the cross sections for 
inelastic and reactive collisions combined. But they only randomize the directed motion 
of the electrons in the electric field. The thereby induced changes of the electron energy 
distribution function affects the plasma-chemistry only indirectly, in as much as an increase 
of energy in the electronic subsystem makes certain collisions more probable than others. 



1.3.1 Production of ions 

The main production processes for ions in a molecular, electro-negative gas discharge are 
electron impact ionization and dissociative electron attachment. Both processes share the 
same compound state, AB~ , but the former leads to positive molecular ions whereas the 
latter to negative atomic ions. Impact ionization is in addition also the main source for 
electrons which, depending on their energy, trigger a multitude of excitation and dissociation 
collisions. 

Electron impact ionization is similar to electronic excitation (see below), except that 
the excited state belongs to the two-electron continuum. In the notation of the previous 
subsection, the differential cross section is thus proportional to the modulus of the reaction 
amplitude squared, 

dai_.p ~ \{^'p^^^jr\¥^^)\'dnMdE2 , (1.35) 

where Qi is the solid angle associated with ki, the momentum of the primary electron, Q2 is 
the solid angle associated with k2, the momentum of the secondary electron, and E2 is the 
energy of the secondary electron. The exit channel state, l^'^k^kj) = |0/3^^^V'kik2)' describes 
the internal state of the positive ion and the six-dimensional relative motion between the 
two electrons and the ion, which is controlled by V, that is, the Coulomb interaction of the 
primary and secondary electron with each other and with the positive ion. The entrance 

channel of the set of Lippmann-Schwinger equations satisfied by the scattering state I^^Qk'') 
contains an inhomogeneity, |$Q,k) = |0a^)|k), where 10^^) is the initial state of the molecule 
and |k) is the plane wave representing the relative motion of the initial electron and the 
molecule. 

From a rigorous mathematical point of view, the exact two-electron scattering state in the 
field of a positive ion is unknown. Hence, exact calculations of electron impact ionization 
cross sections are not available, even for atoms, where complications due to the nuclear 
dynamics are absent. Perturbative treatments, based, for instance, on the distorted Born 
approximation, which replaces "ipk-^^^ by a plane wave for the primary electron and a Coulomb 
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wave for the secondary electron, are possible, provided the velocity of the incident electron 
is much larger than the velocities of the electrons bound in the molecule. Exchange and 
correlation effects can then be ignored and the calculation of cross sections for electron- 
impact ionization reduces essentially to the calculation of photo-ionization cross sections. 
This approach, due originally to Bethe (0) and reviewed by Rudge (25) and Inokuti (27), 
fails at low electron energies. But in that region, Vriens' binary-encounter model (8) can be 
used. It assumes that the primary electron interacts pairwise with target electrons, leaving 
the remaining electrons and the nuclear dynamics undisturbed. The ionization cross section 
is then basically the Mott cross section for the collision of two electrons, appropriately 
modified by the binding and kinetic energies of the target electrons. 

The semi-empirical Kim-Rudd model for electron impact ionization (0) combines the Bethe 
model with the binary-encounter model. For an individual orbital, the ionization cross 
section is then given by 



S 



e+u+1 2 



1 



lne+(2-Q) 



^ 1 In e 

e e + 1 



:i.36) 



with e = E/B and u = U/B, where E is the kinetic energy of the incident electron, U is 
the average kinetic energy of the orbital, and B is the binding energy of the orbital. The 
remaining quantities are S = AnN/B"^ and Q = 2BMf/N with the orbital occupation 
number and Mj an integra l over the differential oscillator strength defined in Ref. (0). The 
results obtained from Eq. (11.361) are surprisingly good. 

Another widely used semi-empirical model for electron impact ionization is the Deutsch- 



Mdrk model ([TO; [ul; UJ) • It uses cross sections for the atomic constituents of the molecules 
and sums them up according to the atomic population of the molecular orbitals which is 
obtained from a Mullikan population analysis. The total ionization cross section can then 
be written as 



(^^(E) = ^gi,nl^^rl^lNi^nlf{e^,nl) 

n,l,i 



[1.37) 



with rf^i the mean square radius of the n, l-th sub-shell of the constituent i, Ni^ni the oc- 
cupancy of that sub-shell, Qi^ni weighting factors which have to be determined empirically 
using reliable atomic cross section data, and a function 



1 / ei,ni - 1 



i,nl) 



3/2 



l + ?(l-^ 

3 \ ICi nl 



In [2.7 + ^ei^ni - Ij 



(1.38) 



which describes the energy dependence of the cross section. The quantity Ci^ni = E/Bi^ni with 
Bi^ni the ionization energy of the considered sub-shell of the i-th constituent. More elaborate 
versions of the Deutsch-Mark model account also for the different angular symmetries of the 
sub-shells (flO| ). 



The cross sections for impact ionization are fairly large because it is an optically allowed 
process (this can be most clearly seen when the Born approximation is applied) and the 
energy in the final state can be distributed among the two electrons in infinite many ways. 
Hence, a gas discharge usually contains a high concentration of electrons and positive ions. 

Negative ions are produced by dissociative attachment, which is a resonant process due 
to the collision of a slow electron with a (in general) vibrationally excited molecule, where 
the compound state AB~ dissociates into a negative ion A~ and a neutral atom B. In 
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Figure 1.6: Schematic representation of the relevant potential energy surfaces for disso- 
ciative attachment. For simplicity, it is assumed that the molecule AB is initially in the 
vibrational groundstate, although this is usually not the case. The electron is thus captured 
at R = Ro, where Rq is the equilibrium distance of the two nuclei (upward directed vertical 
arrow). The thereby created AB~ compound is envisaged to be in an anti-bonding state 
whose potential energy surface crosses at i? = R^ the potential energy surface of the AB 
state. For R < R^ the compound state has thus a finite probability to decay (downward 
directed vertical arrows). However, provided the state survives until the nuclear distance 
R > Rx, it asymptotically reaches the {A~ , B) dissociation limit and dissociative attachment 
is completed. 



principle, it would be also conceivable that the incident electron is permanently bound to 
the molecule and a negative molecular ion is formed. A necessary condition for this to happen 
would be that the relevant potential energy surface of the AB~ state is bonding and has a 
minimum outside the "potential well " corresponding to the electronic groundstate of the 
AB molecule. However, the electron is always captured above the re-detachment threshold, 
that is, the molecular ion is initially in an excited state. It can be only stabilized when it 
looses or disperses excess energy. The former occurs through collisions with other molecules 
or the wall while the latter takes place through vibrational modes of the AB~ compound. 
In most cases, these processes are not very efficient and the AB~ state is only a short-lived 
resonance not affecting the macrophysics of the discharge. 



As pointed out by O'Malley (1131) and others (14), dissociative attachment most likely 
occurs when an anti-bonding state of the collision compound AB~ crosses the initial potential 
energy surface of the AB molecule as schematically shown in Fig. II. 6[ This is the typical 
situation for the resonance model to be applicable. Thus, with the appropriate identifications 
and boundary conditions, t he cr oss section for dissociative attachment can be calculated from 
the effective Hamiltonian (ll.26p . 

The natural boundary conditions for dissociative attachment are an incoming plane wave 
for the relative electron-molecule motion and an outgoing spherical wave in the relative 
(yl~,i?) motion. The formula for the attachment cross section, which is an reaction cross 
section [cp. with Eq. (11.191) ]. would then however contain the interaction in the exit chan- 
nel, that is, the interaction between A~ and B, which is not part of the resonance model. 
This complication can be avoided when the adjoint scattering problem is considered, whose 
boundary conditions are an incoming plane wave in the {A~,B) channel and an outgoing 
spherical wave in the electron-molecule channel. The cross section contains then the interac- 
tion between the resonant state and the e — AB scattering continuum. Specifying Eq. (11.201) 
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to the present case, the differential cross section for dissociative attachment becomes 



da 



DA 



mM K 



(27r)2 



dR 



v,UR)xuAR) 



dQ 



K 



1.39) 



where m and k are the reduced mass and the relative momentum of the (e , AB) system, 
respectively, M and K are the corresponding quantities in the (v4~, B) system, Xv^{R) is the 
vibrational state of the molecule, and '^^^^ {R) is the scattering state satisfying the complex 
conjugate of Eq. fO?!) with Vd{R) = Vab-{R), Vo{R) = Vab{R), Vdk{R) the interaction 
between the anti-bonding AB~ state and the electron-molecule scattering continuum, and 
J{R) = [V^(oo) — V(i{R)] exp[iKR]. The total energy available for the collision E = k'^/2 + 
oju- = K"^ /{2M) + ua- + ojb, where oja- and ojb denote the internal energies of the ion and 
atom, respectively, and ujy. is the vibrational energy of the molecule. 

In order to avoid additional indices, quantum numbers for the internal state of the (A^, B) 
system are suppressed and rotations of the molecule are also ignored. Because the period 
of rotation is much longer than the collision time, rotations could be included within the 
adiabatic approximation. Finally, provided the kinetic energy of the incident electron, A;^/2, 
is much larger than the vibrational energy of the target, uj^-, the local approximation could 
be employed, that is, Eq. (11.341) could be used instead Eq. fll.27p . Further details about 
the calculation of electron attachment cross sections can be found in the review article by 
Chutjian and collaborators (l30l ). 

Bardsley and coworkers ([T3) have shown that within the semi-classical approximation 
da^JdVlK factorizes into a capture cross section, which describes the formation of the com- 
pound state, and a survival probability for that state. The semi-classical calculation more- 
over shows that the capture cross section increases faster with temperature than the survival 
probability decreases. Thus, the attachment cross section increases with temperature. It 
is however always one or two orders of magnitude smaller than the cross section for the 
corresponding elastic electron-molecule scattering. 



1.3.2 Destruction of ions 

Whereas in an electro-negative gas discharge only two processes are primarily responsible for 
the production of ions and electrons - electron impact ionization and dissociative attachment 
- a large number of processes may lead to the destruction of ions and electrons (see Table [LTj) . 
Depending on the parameters of the discharge, the loss of ions may be due primarily to 
recombination or detachment. Recombination may furthermore occur between positive ions 
and electrons or between positive and negative ions (annihilation). For molecular positive 
ions, the former process is usually accompanied by dissociation and is thus called dissociative 
recombination. Ion-ion annihilation, on the other hand, is in most cases only a charge 
transfer, which does not lead to a rearrangement of the nuclei. Finally, detachment of 
negative ions can be either initiated by electrons or by neutrals. Electron-induced detachment 
is similar to electron impact ionization whereas detachment due to (molecular) neutrals may 
lead to dissociation, as well as, when the neutral is in a meta-stable state, to association. 

Dissociative recombination is triggered by slow electrons. It is thus a resonant process, 
similar to dissociative attachment. Very often, it is the dominant loss process for positive 
ions and electrons with relatively large cross sections because at least one of the atomic 
fragments in the exit channel is usually in an excited state implying that the energy in the 
exit channel can be distributed in many different ways (|3^ . 
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Figure 1.7: Dissociative recombination occurs when an anti-bonding potential energy sur- 
face of the (excited) molecule, Vab{R), crosses the potential energy surface of the positive 
molecular ion AB^ . Although the states involved are different, the mechanism is similar to 
dissociative attachment. When the molecule, which is the collision compound in this case, 
survives auto-detachment for R < R^, it may reach the {A, B) dissociation limit. As a result, 
dissociative recombination takes place. 



The resonance model for dissociative recombination, originally proposed by Bardsley (jl 



18l ). is based on the potential energy surface diagram shown in Fig. I1.7[ The most favorable 
situation for the process to take place is when an anti-bonding potential energy surface 
of the molecule AB., which is the collision compound for this reaction, crosses the potential 
energy surface of the positive ion and is for large inter-nuclear distances below the vibrational 
ground state of the ion. This is, for instance, the case for H2 , N2 , and O^. 

The cross section for dissociative recombination can then be casted into the form fll.39p . 
with M and K the reduced mass and the relative momentum of the {A, B) system, respec- 
tively, m and k the corresponding quantities of the {e~,AB~^) system, Xui{R) the initial 

vibrational state of AB~^, and "^[^{R) the solution of the complex conjugate to Eq. (11.271) 
with Vd{R) = Vab{R), Vo{R) = Vab+{R), and Vdk{R) the interaction between the AB state 
and the — AB^ scattering continuum. The inhomogeneity representing the boundary 
condition, J{R) = [Vd{oo) — V^(i?)]$^(i?), contains now a Coulomb wave, $^(-R), because 
the scattering continuum is for two charged particles: an electron and a positive ion. More 



sophisticated approaches use quantum-defects to characterize the scattering states ( I43l : |4 

In electro-negative gas discharges, there is an additional recombination channel: Ion- 
ion an nihil ation. The potential energy surface diagram relevant for this process is shown 
in Fig. 11.81 The essential point is that at an inter-nuclear spacing R = R^ the potential 
energy surface for the {A~, AB~^) configuration, which decreases with decreasing R because 
of the Coulomb attraction between the two ions, crosses an anti-bonding potential energy 
surface of the A2B collision compound. At this separation, the probability for the system 
to switch from the {A~,AB~^) to the {A,AB) configuration is particularly high. Based on 
this observation, a semi-empirical Landau-Zener model can be constructed which relates the 
cross section for ion-ion annihilation to the ionization energy of the AB molecule and the 
electron affinity of the A atom (see Fig. II. 8p . 

The Landau-Zener model illustrates quite nicely how semi-empirical models encode com- 
plicated processes in a few physically intuitive parameters. In contrast to the processes 
discussed so far, where the relevant nuclear dynamics took place on a single potential energy 
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Figure 1.8: Relation of the potential energy surfaces typical for ion-ion annihilation. In the 
entrance channel, the nuclear dynamics is governed by l^(-R), the potential energy surface 
of the {^AB^^A~^ configuration, whereas the dynamics in the exit channel is controlled by 
V/(-R), which is the potential energy surface of the (Ai?, A) configuration. The strong mixing 
of the two configurations encoded in V[K) (see Eq. f ll.42p ) is maximal at R = R^- Both, 
Vi{R) and Vf{R), belong to the collision compound A2B. The potential energy surface for 
the ^2-8"*" molecule is also shown. It allows to define the relevant energies: the ionization 
energy of the AB molecule, E:^, and electron affinity of the A atom, E^. 



surface, ion-ion annihilation forces the nuclei to switch between two potential energy surfaces 
of the collision compound. The minimal theoretical model is therefore based on two coupled 
Lippmann-Schwinger equations for the relative motion of the nuclei. The Landau-Zener 
model is the semi-classical approximation to this set of equations. Following Olson (t20l ). the 
total ion-ion annihilation cross section is then given by 



a"^(E) = AnRi 



AE' 
~E~ 



E{X) 



with 



and 



F(A) 



dx exp(— Ax) 



1 — exp(— Ax) 



A = 27r 



M 




\V{R.) 



2 \V^ -ViWE + AE 



;i.4o) 



;i.4i) 



(1.42) 



where E and M are the kinetic energy and the reduced mass of the relative motion of the 
{A~,AB~^) system, AE is the energy gain due to annihilation, V{Rx) is the interaction 
between the {A~,AB~^) and {A,AB) configurations at i? = Rx, and V/j = dVij{Rx)/dR 

with Vi{R) and Vf{R) the potential energy surfaces of the {A~,AB'^) and {A,AB) system, 
respectively. 

Equation (11.401) can be developed further, by recalling that Vi{R) ~ R^^ (Coulomb inter- 
action between A~ and AB~^) and Vf{R) ~ with n > 1 (polarization interaction between 
A and AB). Hence, for large enough Rx, \V- - Vj] ^ R~'^ and AE ^ -Ui{Rx) = R~^. Usu- 
ally, AE > E. Combining all this leads to A ^ ^/2MtiRI'^\V{Rx)\'^ (HO). Since F(A) 
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approaches its maximal value Fmax ~ 0.1 at Xmax ~ 0.424, a rough estimate for the annihi- 
lation cross section is 



1 

1 + 



RxE 



;i.43) 



with Rx determined from Xmax = V2M7rRl^^\V{Ra:)\^ or from empirical cross section data for 
high energies where E ^ R~^ and — > 1.3R^. To determine V{Rx) is not trivial. Ideally, 
it can be parameterized in terms of an effective ionization energy Ef^^ of the AB molecule and 



the electron affinity E^ of the A atom. Olson and coworker (|21| ) obtain for large inter-nuclear 
distances V{R) = 1.0U^y Ef^E^^R exp[-0. 857 R] with R = 0.5{^/2E^ + ^2Ef^)R and 
q the Franck-Condon factor that represents the overlap of the vibrational states (see below). 
In order to apply this formula, the ionization potential and the electron affinity have to be 
known. For atomic negative ions this is no problem. But for molecular negative ions this 
can be subtle because the ion may be in an (unknown) excited state. Fortunately, molecular 
negative ions have a rather short lifetime, as discussed before. They are thus not of major 
concern in the commonly used gas discharges. 

For molecular positive ions, dissociative ion-ion annihilation may be also possible. It 
contains two relative motions in the exit channel and is thus much harder to analyze the- 
oretically. Experimentally, it is not always possible to discriminate between this process 
and the charge-transfer-type reaction discussed in the previous paragraphs. The energy de- 
pendence of the two cross sections is however essentially the same, except that they reach 
the constant value at different energies. Thus, when R^ is determined from the high en- 
ergy asymptotics of empirical data, the annihilation cross section is in fact an effective cross 
section, comprising both charge-transfer-type and dissociative ion-ion annihilation. 

In contrast to positive ions, negative ions can be also destroyed by slow collisions with 
neutral particles (detachment). To be specific, associative detachment oi negative ions due 
to molecules is discussed. Microscopically, it is a resonant process similar to dissociative 
attachment and dissoci ativ e recombination. An auto-detaching, resonant state plays thus an 
important role. In Fig. II. 91 the relevant potential energy surfaces are shown. The asymptoti- 
cally stable configuration (y4~, AB) changes with decreasing R into a quasi-bound, resonant 
collision compound, A2B~ , from which the electron auto-detaches, leaving behind a free 



electron and a neutral particle in an excited vibrational state (lla : \Va \m \ 



The differential cross section for this process can be obtained from the resonance model, 
using as boundary conditions an incoming plane wave in the {A~,AB) channel and an 
outgoing spherical wave in th e (e~, A2B) channel. Normalizing continuum states on the 
energy scale and applying Eq. (11.191) leads to 



da 



AD 



d^k^. , (1.44) 



where Xuf{R) is the vibrational state of the A2B molecule and "^^^ satisfies Eq. (11.271) 
with Vd{R) = Va2B-{R), Vo(-R) = K42b(-R), and Vdk{R) the interaction between the quasi- 
bound ^2-8" state and the e~ — A2B scattering continuum. The inhomogeneity is J{R) = 
[Vd(oo) — Vd{R)]^K{R) with ^k{R) a plane wave for the relative motion of the {A~,AB) 
system. Energy conservation enforces E = K'^/ {2M) + Vd{oo) = k"^ / {2m) + uj^j where M and 
m are the reduced masses of the {A~.,AB) and the {e~,A2B) system, respectively. Notice, 
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Figure 1.9: Illustration of associative detachment. The potential energy surface of a bond- 
ing state of the collision compound ^2-8" crosses the potential energy surface of the A2B 
molecule. Auto-detachment of the compound state produces then an electron and a vibra- 
tionally excited ^42-8 molecule. Also shown is an anti-bonding A2B~ state (dashed line). 
Detachment may also occur from such a state, but then the initial kinetic energy of the 
{A~ , AB) system has to be larger than Va2B-{.R'x) ~ K42B-(^)- 



since da^J^ is a differential cross section for reactive scattering, the interaction in the exit 
channel, Vdk, appears in Eq. fll.44p . as it should be. 

The most favorable situation for associative detachment is when the potential energy 
surface of the neutral molecule in the exit channel supports a bound state whose dissociation 
energyE'jf is larger than the electron affinity E-^ of the A atom. This situation is shown in 
Fig. 11.91 Provided detachment is mediated by an attractive state of the collision compound, 
it takes place even for vanishing initial kinetic energy in the {A~ , AB) channel. If this is the 
case, detachment is a very efficient loss channel for negative ions, even at low temperatures. 
If, on the other hand, the compound state is anti-bonding, associative detachment occurs 
only when the initial kinetic energy of the colliding particles is larger than Vd{R'x) — Vd(oo), 
where R!^ is the point where the repulsive potential energy surface crosses Vo(-R). 



1.3.3 Excitation of internal degrees 



The chemistry and charge balance of a gas discharge is also affected by inelastic collisions, 
that is, collisions which increase the internal energy of molecules, atoms, and ions. Excited 
(meta-stable) particles are reactive and participate in basically all particle number changing 
collisions. In low temperature gas discharges, vibrationally excited molecules play a particu- 
larly important role, because, for typical operating conditions, they are efficiently produced 
by resonant electron-molecule scattering. 

The kinetic energy of electrons in a low-temperature gas discharge is typically a few 
electron volts. At these energies, the electron-molecule collision time is rather long, favoring 
therefore resonant enhancement of the collision. The cross section for vibrational excitation 
of molecules is thus given by Eq. (11.141) with the scattering state obtained from the resonance 
model. Specifically for continuum states normalized on the energy scale, the cross section 
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Figure 1.10: Illustration of resonant vibrational excitation of a molecule for a large (a) and a 
small (b) imaginary part of F{R, R'; E). The dashed line indicates the dissociation threshold, 
Vo(-R) is the potential energy surface for the molecule, Vd{R) is the potential energy surface 
for the auto-detaching resonant state of the negatively charged molecule, and the insets show 
the typical shape of the cross section for the two cases. 



becomes 

2 

dVlf , (1.45) 

where Xvf{R) is the vibrational state of the molecule after the collision and '^^^^{R) is 
the solution of the effective Lippmann-Schwinger equation fll.27p with Vd{R) = Vab-{R), 
Vo(-R) = Vab{R), Vdk{R) the interaction between the resonant AB'^ state and the electron- 
molecule scattering continuum, and J{R) = Vfiin{R)Xui{R)- Energy conservation implies 
E = kf/{2^) + ujy. = k'j/{2fi) + ojy^ with /i the reduced electron-molecule mass. 

The shape of the vibrational excitation cr oss se ction depends on the imaginary part of 
F{R, R'; E), that is, on T{R, R'; E) (see Eqs. ([L2H]) and fOOD l When it is large compared 
to the inverse of the period of vibration in the resonant state, that is, for a short lived 
resonance, the cross section is smooth. For a long-lived resonance, that is, when T{R, R'; E) 
is much smaller than the period of vibration, the cross sections consist of a series of peaks 
reflecting the vibrational states of the resonant stat e. Th is dependence of the cross section 
on the lifetime of the resonance is illustrated in Fig. ll.lOi 

When the energy exceeds the dissociation energy, dissociation may occur, either indirectly 
through the resonant state, or directly through a transition to the continuum of the nuclear 
motion of the molecule. The indirect process contributes only when the target is initially 
in a high excited vibrational state (l45f). Usually, the direct process dominates. For large 
enough energies, the adiabatic approximation applies and the dissociation cross section is 
given by 

= ^^\T^^.^X^f,Q,)\'dQf (1.46) 



dRxUR)v:^(R)^^;XR) 
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with 

T^.^.iQf, Qi) = J dRF:^{R)t{nf, Qi, R)F,,{R) , (1.47) 

where t{Qf,Qi,R) is the electronically elastic fixed-nuclei scattering amplitude and Fi^^.{R) 
is a continuum nuclear wavefunction of the molecule. Electronic quantum numbers are 
suppressed because they are unchanged; vibrational excitation and dissociation involve only 
one potential energy surface of the target molecule. 

Finally, the production of electronically excited molecules occurs due to scattering with 
high energy electrons. For energies far away from the dissociation threshold, the adiabatic 
approximation holds. Therefore, ignoring the rotational degrees of freedom, the transition 
amplitude can be written as 

Tnf^fniiyii^f,^i) = J dRF*^^^{R)tnfni{^f,^i,R)FniUiiR) 

= tn,nA^f,ni,Ro) jdRF:^,^{R)F^,,,{R) . (1.48) 

In the second line the fact has been utilized that Fn^„- (R) is strongly peaked at R = Rq, where 
Rq is the position of the minimum of the potential energy surface of the initial electronic 
state of the molecule. The differential cross section for electronic excitation is thus given by 



da: 



da3^_ 



EE _ ni-*nf 



f 



J dRF:^,^{R)F^,,XR) 



2 



dQf , (1.49) 



where the first factor is the fixed-nuclei differential cross section &X, R = Rq for the electronic 
transition rii —>■ rif and the second is the Pranck-Condon factor q already mentioned in 
connection with ion-ion annihilation. 



1.4 Concluding remarks 

This chapter discussed elementary collision processes, as they typically occur in molecular 
gas discharges, focusing, in particular, on inelastic and reactive collisions which change 
the internal energy and composition of the scattering fragments: ionization, attachment, 
recombination, annihilation, detachment, and excitation. These processes are the primary 
driving force of the plasma-chemistry in these discharges. Elastic scattering, on the other 
hand, is not explicitly discussed, because its main effect, depositing energy into the electronic 
subsystem, affects the plasma-chemistry only indirectly. 

Instead of merely listing kitchen-made cross section formulae and unrelated cross section 
data, emphasis has been given on an unified description of elementary processes based on 
general principles of quantum-mechanical multi-channel scattering theory. By necessity, the 
presentation is rather dense. Technical details left out, as well as cross section data, can be 
found, respectively, in the original papers and the review articles, monographs, and web-sites 
mentioned in the introduction. 

As far as the kinetic modeling of gas discharges is concerned, collision cross sections 
are input parameters for the collision integrals of the Boltzmann equations. They are thus 
always convoluted with distribution functions. Some of the quantum-mechanical intricacies 
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of particle en- and de-tanglement occurring during a collision are thus eventually averaged 
out. Calculating cross sections earmarked for plasma-chemical applications with the most 
sophisticated quantum-mechanical methods is thus not only not feasible for the (poly-atomic) 
gases of interest, it would be also overkill. Cross sections obtained from effective (semi- 
empirical) models containing only those microscopic degrees of freedom which, for given 
external control parameters, may eventually become active in the collision integrals, should 
be in fact sufficient. A systematic effort to develop (and solve) this type of models for the 
various processes may significantly increase the predictability of plasma modeling. It may 
thus help to turn plasma processing as a whole from an art to a science (46). 
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